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Influence of Spatial Correlations in Strongly Correlated Electron Systems: 
Extension to Dynamical Mean Field Approximation 

Hiroaki Kusunose* 
Department of Physics, Tohoku University, Sendai 980-8578 

We propose a formalism to take account of the correction of the spatial fluctuations to the local self- 
energy obtained by the dynamical mean-field approximation. For this purpose, the approximate dynamical 
susceptibility in the framework of the iterated perturbation theory is proposed and examined. Using the 
formalism, it is demonstrated that the one-particle spectral intensity in the two-dimensional Hubbard model 
at half-filling exhibits the pseudo-gap behavior in the central coherent quasiparticle peak due to the critical 
antiferromagnetic fluctuation. The specific heat is considerably enhanced by the short-range order, which 
assists a tendency of the Mott localization showing the reduction of the double occupancy. We briefly 
discuss a formulation for the superconducting transition temperature in the present approximation. 

KEYWORDS: dynamical mean field approximation, spin fluctuation, metal-insulator transition, Hubbard 
model, iterated perturbation theory anisotropic superconductivity 



1. Introduction 

Dynamics of strongly correlated electron systems has 
been investigated extensively in context of low-temperature 
physics in the heavy fermion compounds, the transition met- 
als and some organic materials. One of the key feature is a 
duality of electronic states in one-particle density of states 
(DOS). Namely, the high-energy incoherent states (so-called 
the Hubbard bands) have a spatially localized character, while 
the low-energy states exhibit an itinerant character with strong 
temperature dependence. Particle-hole excitations within the 
upper or lower Hubbard bands constitute the localized spin 
degrees of freedom, which govern the low-energy spin dy- 
namics. This low-energy spin fluctuation in the two-particle 
correlation can have a significant influence on the low-energy 
itinerant quasiparticle state.'"'' 

Various theoretical approaches have been devoted to inves- 
tigate such an interesting but complicated dynamics. The rep- 
resentative theories among them may be classified by their 
nature of the irreducible vertices, i.e., a source of the ef- 
fective two-body interactions. They depend on the incom- 
ing and the outgoing momenta and energies in general.'* The 
random phase approximation (RPA) approximates the irre- 
ducible vertices to bare interactions, while one-particle states 
are considered in the mean-field level. Reflecting the RPA- 
type two-particle fluctuations to the one-particle self-energy 
in a self-consistent manner, we obtain the fluctuation ex- 
change (FLEX), or the Migdal-Eliashberg approximation.^-^ 
This is a conserving approximation in the sense of Baym- 
KadanofF,'"^ but it violates a number of crucial sum-rules.'" 
On the contrary, the two-particle self-consistent (TPSC) the- 
ory respects primarily the two-particle fluctuations and their 
sum rules rather than the one-particle quantities.'" Following 
this philosophy, local and instantaneous (contact-type) irre- 
ducible vertices are introduced independently in the charge 
and the spin channels, which are determined so as to satisfy 
the local-moment sum-rules in each channel. The one-particle 
self-energy is then improved by including the effect of the col- 
lective modes. All of the approaches mentioned above is the 
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weak-coupling theory in the sense that they neglect energy 
dependences in the UTeducible vertices. Therefore, they fail 
to describe the spin dynamics of the localized spin degrees of 
freedom in strongly correlated electron systems. The renor- 
malized energy scales such as a superexchange interaction 
is never involved in the two-particle correlation. Neverthe- 
less, it is worth noticing that the weak-coupling theories using 
the contact-type irreducible vertices qualitatively account for 
the characteristic momentum dependences inherent from the 
electronic band structures." 

On the other hand, the dynamical mean-field theory 
(DMFT) has been proven successfully to describe the dual 
nature of the electronic states in the energy space, e.g., the 
opening of the Mott gap at half-filling for sufficiently strong 
Coulomb repulsion.'^' '"^ If we remind the fact that the self- 
energy is expressed in terms of the vertices,''* the proper en- 
ergy dependences in the irreducible vertices have crucial im- 
portance in describing the duality of the one-particle elec- 
tronic states as well as the renormalized energy scale in 
the two-particle correlations. The one-particle self-energy in 
DMFT does not contain any contributions from the nonlocal 
coiTelations. Hence, extensive studies such as cluster exten- 
sions of DMFT have been carried out to take account of short- 
range correlations.'^"'^ 

In this paper, we propose an alternative extension to DMFT 
by taking account of the coiTections of the spatial fluctuations 
in the self-energy. Sadovskii and co-workers have proposed 
an introduction of the effective length scale of the param- 
agnon excitations for this purpose.''' -" We assume local but 
energy-dependent irreducible vertices to construct the nonlo- 
cal two-particle correlations. This is a minimum requirement 
to describe the two-particle coiTelations both with renormal- 
ized energy scale and band-specific features in the momen- 
tum space. The resultant one-particle states with the nonlocal 
self-energy can represent the duality in the energy space with 
a proper momentum dependence. A similar approach named 
as the dynamical vertex approximation (DFA) has recently 
been developed on the basis of the parquet equation.^' Specif- 
ically, we extract the local irreducible vertices from the re- 
sult of DMFT, and construct the nonlocal dynamical suscep- 
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tibilities. Then the nonlocal corrections of the spatial correla- 
tions are taken into account to improve the DMFT self-energy 
analogous to the TPSC theory. For this purpose, we develop 
an approximate construction for the local dynamical suscep- 
tibiUty based on the iterated perturbation theory (IPT).^^^^^ 
Then, we discuss the influence of the spatial correlations in 
the one-particle spectral intensity for the symmetric 2D Hub- 
bard model at half-filling. 

The paper is organized as follows. In §2, we estabUsh the 
formulation of the extension to DMFT. The subsequent sec- 
tion is devoted to explain the construction of the dynamical 
spin susceptibiUty based on IPT. The validity of the approx- 
imate spin susceptibihty and the U-T phase diagram are dis- 
cussed. The main result of the paper is given in §4. The influ- 
ence of the spatial correlations in the one-particle spectral in- 
tensity is mainly discussed. The summary and discussion are 
presented in the last section. The outiine of the IPT scheme 
extensively used in this paper is summarized in Appendix A. 
A formulation for the superconducting transition temperature 
based on the present approximation is briefly discussed in Ap- 
pendix B. 

2. Extension to DMFT 

In this section, we describe the formulation to take into ac- 
count the correction of the spatial correlations to the DMFT 
self-energy. To be specific, we consider the Hubbard model 
on the two-dimensional square lattice. 



(2.1) 



where cj[^ is a creation operator of conduction electron with 
the wave- vector k = (kx, ky) and the spin cr. The electron dis- 
persion is given by et = -2t{coskxa + cos kyd). The overlap 
integral t and the lattice constant a are set by / = a = 1. 
tiia- = cJ^Cia- is the density operator at the site i. For no- 
tational simplicity, we use the (2+l)-vector representation, 
k = (k, icOn), if necessary. The 3-vector summation, J]^, repre- 
sents {TIN) Y,k 2n' ^ being the number of sites. We also use 
the notational conventions, (i) the letters k and w„ {q and e„) 
are used to represent the fermionic (bosonic) Matsubara fre- 
quency, and (ii) a quantity computed by DMFT is indicated 
by the superscript, L. We consider only the symmetric case 
with the electron density, n = 1, throughout this paper. 

We begin with two exact relations. One is the Bethe- 
Salpeter (BS) equation in the charge and the spin chan- 
nels,4-8-io 



-Yj^cbik,k";q)xa.ik",k';q) 

k" 



Xsp(k,k';q) ^xo(k;q) 



+ Y,^svik,k";q)xspik",k';q) 



(2.2a) 



(2.2b) 



where xo(k', q) - -G(k)G(k + q) is the lowest-order two- 
particle Green's function. Xchik, k'; q) and Tcbik, kf; q) are the 
two-particle Green's function and the irreducible vertex in the 
charge channel. The similar notation is applied for the longi- 
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Fig. 1. The Feynman diagrams for (a) the Bethe-Salpeter equation, (2.2), 
(b) the self-energy in terms of susceptibilities and irreducible vertices, 
(2.5). 



tudinal spin channel. The sign of the irreducible vertices are 
chosen so as to become the positive U in the lowest order, i.e. 
rch(fc, k';q)^ rsp(yt, k';q)^U for RPA or FLEX approxima- 
tions. Note that rch(^, k';q) - Uch and Fspik, k';q) - Usp are 
used in TPSC. The ordinary susceptibility is obtained by the 
summations, 



Xa(q) = ^ Ara(^. ?). (« = ch, sp). (2.3) 

k,k' 

Note that we have defined the (spin) density operator as 
Pch(sp) - (c^^c-i+c^j^ci)/ in the susceptibilities. TheFeyrmian 
diagram for the BS equation is shown in Fig. 1(a). 

The Dyson's equation for the one-particle Green's function 
is given by 



G-\k) = ioJn+ri-ek-^{k), 



(2.4) 



where Y.{k) is the self-energy in which the Hartree contribu- 
tion is subtracted. Namely, the chemical potential, fi, con- 
tains the Hartree contribution as // = yu - f/n/2 where n = 
2 Tik G{k)e'""^* . We omit spin indices since we will discuss 
only the paramagnetic state. The self-energy can be expressed 
in terms of the charge and the (longitudinal) spin susceptibil- 



Uk) ^^Yj Z[^^''(^, k"; q)xc^ik", k'-q) 



q k',k' 

+ r,p(fe, r ; q)Xsp(k", k'; q)\Gik + q). (2.5) 

This is the second exact relation, and its Feynman diagram 
is shown in Fig. 1(b). It should be emphasized that we can 
also write down the self-energy in terms of alternative fluctu- 
ation, such as superconducting (pairing) one.^^ Here, we have 
chosen the expression of the charge and the spin fluctuations, 
since the latter fluctuation plays a dominant role in strongly 
repulsive U regions. For negative U problem, the expression 
of the pairing fluctuation is appropriate.'^^ 

In order to solve the above two equations, let us intro- 
duce two approximations as follows. As was mentioned in 
the introduction, the energy dependence in the irreducible ver- 
tices is significantly important in strongly correlated electron 
systems. Therefore, as the simplest approximation, we retain 
only the external energy dependence, in the irreducible 
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vertices, 

r„ik,k';q)~TM- (2.6) 
Then, the BS equations, (2.2), can be solved formally as 

x;^iq)=Xo\q) + rcUem), (2.7a) 
x;p(q)=Xo'(q)-^sp(€„). (2.7b) 

Note that the momentum dependence is considered only 

through ;yfo(<7). Within the same approximation, we can divide 
the contributions from the fluctuations in eq. (2.5) into the lo- 
cal and the nonlocal ones, 

U 



tik) = l{oJn) + - 2[rch(em)A;trch(^) 



+ rsp(e„,)A^sp(9) 



G(k + q), (2.8) 



where t^Xaiq) = XaiM) -Xai^m) describes the nonlocal contri- 
bution. The local self-energy is defined as 



S(w„) 



+ rsp(e„);kfsp(em) G((o„ + e^,), (2.9) 
where the local Green's function is given by G(a>„) = 

Next, we introduce the second approximation. The local 
quantities in eqs. (2.7) and (2.8) are replaced with those ob- 
tained by DMFT, 

£(w„) ~ tHco„), 



Xai^m) ~ Xa^^in^t 



(2.10) 
(2.11) 



and the local vertices are calculated from the relations similar 
to eqs. (2.7), 

+ reh(em) ~ 1 Uchi^m) - 1 /X^^iCm), (2. 12a) 

- rsp(e„) ~ 1 /x%(em) - 1 lx\Mm\ (2. 12b) 

where the local irreducible susceptibility is given byxo(em) = 
-T 2(j„ G^{ojn)G^{cOn + em)- With these approximations, (2.8) 
reads the equation to improve the local DMFT self-energy by 
including the corrections of the spatial correlations. Then, the 
updated Green's function is obtained by the Dyson's equation 
with the improved self-energy. The chemical potential fx is 
determined such that 2 Y^k G(^)e'""°+ coincides with the given 
electron density.'" In this paper, we always use jx = Q because 
of the symmetric condition. 

3. Dynamical Susceptibility Based on IPT 

In the previous section, we propose the formulation to im- 
prove the local self-energy by including the correction of the 
nonlocal bosonic fluctuations. To perform a practical calcula- 
tion, the dynamical susceptibilities as well as the local self- 
energy are required within DMFT. In practice, however, effi- 
cient algorithms to compute the dynamical susceptibilities at 
finite temperature are very limited. In fact the quantum monte 
carlo (QMC) simulation based on the Hirsch-Fye algorithm 
is only such an efficient algorithm at moment.''^'^^''^' In this 
section, we first develop an approximate construction of the 
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dynamical susceptibility based on IPT. Then, we examine the 
vahdity of the approximate dynamical susceptibility. 

3.1 Formulation of dynamical susceptibilities 

Within the spirit of DMFT, we have the local version of the 
relations, 2^'^° eqs. (2.2) and (2.5), 



(3.13a) 



xLioJn, oj'„; em) = x^ioj„; e„) 



(3.13b) 



and 



tH<On) = ^T'YjTj [rch('^-<; e„)x^^ico':,aj'„; e.) 



-I- r^p(w„, w"; e,„)x^p(w;', w;; e,„) 



G^o;„ + e,„). (3.14) 



Eliminating the irreducible vertices from these relations, we 
obtain 



^Wn) = -^T Yj[x^p(^n; em) - ^t^bi^n, ^m) 



X GHoJn + em)/xo(oJn; e„,), (3.15) 

where ^^(w„; em) = T Y,co'„X^r(<^n,(^'„;em). This expression 
implies that the local self-energy contains information of the 
two-particle susceptibility (it is nothing but the spin antisym- 
metrized contribution). 

Meanwhile, the self-energy in IPT (see Appendix A for de- 
tail) has the form, 

^H<On) = 1^1] K(cj„)Mem)0((^n + €m), (3.16) 

with K{a>n) = 2UAI[\ - Bl!^H(^„)] and xd^,,,) = 
-T Yjoj,, ©{(^n)@{oin + em), where ll^H^con) being the second- 
order self-energy with respect to the cavity Green's function, 
QitOn)- Comparing this expression with eq. (3.15), we obtain 

4)(fm) - A4(fm) = -Mem)T ^ K{cOn)G\lOn) 

CJfi 

X [g(oj„ + e,„) + §(0Jn - e„,)]/2. (3.17) 

Note that;^^^p(e„,) is written in the symmetrized form with re- 
spect toe,,, «-> -e,„. This expression can be calculatedby using 
only the one-particle quantities in IPT. 

In order to separate the charge and the spin contributions, 
we further approximate r^j^(em) ~ t/ch, which will be deter- 
mined so as to satisfy the sum-rule,'" 

T ^A^^^-) = \ - «) + ' (3.18) 

with Xch(^'») = Xo(fm)/[l + UcbX^iem)]. The double occu- 
pancy, d = (ni^nn), is computed by IPT, i.e., d = (n/2)^ -i- 
(T/U) 2:„„ tHco„)GHoJ„). Then, using eq. (3.17) andAr^^(em) 
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Fig. 2. The T dependence of the static local spin susceptibility for U = 
10 and 16. The dotted line represents the Curie law in the atomic Umit. 




Fig. 3. The U dependence of the static local spin susceptibility for T = 0.1 
and 0.2. The arrows indicate the hysteresis loop. 



with t/ch, we can compute the approximate local spin suscep- 
tibility within the IPT scheme. The charge vertex Ucb is turned 
out to be rapidly increase function of U, and hence the charge 
susceptibility is negligible in strong U region. 

3.2 Validity of approximate spin susceptibility 

Let us examine the validity of the approximate spin suscep- 
tibility as follows. We have used the tight-binding DOS, po(e), 
of the two-dimensional square lattice with the 2048 Matsub- 
ara frequencies in solving the DMFT equations. For the sum- 
mations of the convolution type, the use of the Fast Fourier 
Transform (FFT) algorithm is quite efficient.'''^' Retarded 
quantities are obtained by the numerical analytic continuation 
to the upper half plane, i.e., iun — > co + i6{ie„ — > e+iS), using 
the Pade approximant.^^-^^ 

The T dependence of the static local spin susceptibility for 
f/ = 8, 10 and 16 is shown in Fig. 2. The dotted line repre- 
sents the Curie law in the atomic limit, {x\^{€m = 0) = 1 12T). 
The T dependence of ;if^p(0) approaches the Curie law as U 
increases. It is important to note that a simple irreducible ver- 
tex based on a weak-coupling perturbation theory cannot re- 
produce the behaviors in the atomic limit (the relevant energy 
scale, /, is inversely proportional to U). The discontinuous 
jump at r ~ 0.18 for U = 10 is due to the metal-insulator 
transition (see the U-T phase diagram in Fig. 7). Such a sin- 
gular behavior is confirmed by the QMC calculation.^^ 

The U dependence of x]^p(0) for T - 0.1 and 0.2 is shown 
in Fig. 3. ;if^p(0) for T = 0.2 rapidly increases around the 
metal-insulator crossover region. Two discontinuous jumps 
for r = 0.1 are due to the first-order metal-insulator transi- 
tion. The hysteresis loop of the first-order transition is indi- 
cated by the arrows. These behaviors are also confirmed by 
the QMC calculation. '2 

The spectral intensity ofx^p(e) for f/ = 8 is shown in Fig. 4. 
The low-energy weight of the spectral intensity develops as T 
decreases. Notice that the peak of the intensity is located at 
the order of the superexchange interaction, J ~ t^/U, irre- 
spective of the presence of the quasiparticle coherent states at 
CO = (see, Figs. 7 and 9b). This observation supports the as- 
sertion emphasized by Vilk and Tremblay'" that the bosonic 
collective modes are weakly dependent on the precise form of 
the single-particle excitations. On the other hand, the details 
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Fig. 4. The spectral intensity of the local spin susceptibility for f/ = 8 and 
T = 0.6, 0.4, 0.3 and 0.25. 

of the single-particle self-energy can be strongly influenced 
by scattering from collective (paramagnon) modes. Hence, 
we first compute the two-particle correlations with the DMFT 
Green's functions, and then we improve on the self-energy by 
including the effect of collective modes as shown later. 

Using eqs. (2.12) and (2.7b), the momentum dependence of 
the spin susceptibihty is evaluated. Figure 5 shows the spec- 
tral intensity of Xsp(Q,e) for f/ = 8 at Q = {n,n). Because 
of the perfect nesting of the Fermi surface, the low-energy 
weight of the staggered component is strongly enhanced as 
T approaches the antiferromagnetic (AF) instability, (see 
Fig. 7). The inset shows the average of Im;^^sp(^', e) over the q 
space. As compared with Fig. 4, the slope of the local spectral 
intensity at e = becomes much steeper toward the AF insta- 
bility due to the critical divergence of the staggered spin sus- 
ceptibility. It results in the strong enhancement of the NMR 
relaxation rate, l/TiT = lim^-^o lmxspie)/e. This critical spin 
fluctuation simultaneously causes the pseudo-gap behavior in 
the one-particle DOS as will be discussed later. Thus, the 
pseudo-gap behavior due to the spin fluctuation should always 
be accompanied by the enhancement of \ I T\T rather than a 
spin-gap behavior. 

The contour plot of the spectral intensity of ;^sp(?. along 



J. Phys. Soc. Jpn. 



Full Paper 



Hiroaki Kusunose 5 



10 



X 

S 5 



T=().25 




T=0.3 



T=0.4 



U=8 



T=0.6 



Fig. 5. The spectral intensity of the spin susceptibility at Q = {7t,7t) for 
{/ = 8 and T = 0.6, 0.4, 0.3 and 0.25. The inset shows the spectral intensity 
of the local spin susceptibility. 



U=8. T=0.3 




Fig. 6. The contour plot of the spectral intensity of the spin susceptibility 
along the high-symmetry lines for U = S and 7" = 0.3. 



the high-symmetry Unes for f/ = 8 and T = 0.3 is shown in 
Fig. 6. The high-symmetry points are defined as F = (0, 0), 
X = (7r,0), M = (n,7T) and M' = (n/2,n/2). It exhibits the 
precursor of the collective spin-wave-like dispersion of the 
energy-scale J with the highest intensity at Q. The spin fluc- 
tuation causes a significant modification on the DMFT self- 
energy, as will be shown in the next section. 

Now, we evaluate the Neel temperature, by the condi- 
tion, XspiQ' 0) - 10"^. Note that due to the approximate spin 
vertex, the critical behavior has the mean-field character with 
finite Tn even in two dimension. The U-T phase diagram is 
shown in Fig. 7. Td and Td are the first-order metal-insulator 
transition lines determined by IPT. We define the coherent 
temperature, Tcoh, at which the central quasiparticle peak be- 
gins to appear in the DMFT DOS, A(w) = -lmG^(Lo). With 
increase of U, Tn first increases, showing the maximum at 
U ~ 7, and then decreases followed by 2/U. The suppression 
of Tn around U - 10 is ascribed to the singular behavior of 
the static local susceptibility, ;t^^p(0), due to the metal-insulator 
transition (see Figs. 2 and 3). In other words, the drop of 
x]^p(Q) leads to the suppression of the local spin vertex via 
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Fig. 7. The U-T phase diagram for half-filling. T^i and Ti.2 are the metal- 
insulator transition temperatures determined by DMFT. r^oi, is the coher- 
ent temperature, below which the coherent state develops near the chemical 
potential. The Neel temperature Tn is determined by eq. (2.7b). 



U=8, T=0.245 




Fig. 8. The contour plot of the spectral intensity A(k, to) along the high- 
symmetry lines of the Brillouin zone for f/ = 8 and T = 0.245. The 
pseudo-gap behavior is observed at the X-M' Une. 



eq. (2. 12). In comparison with the result of the QMC'^'^^ with 
the full local vertex, r^p(aj„, w^; e,„), the present Tn is over- 
estimated in weak-coupling region, and underestimated in 
strong-coupling region. The source of the discrepancy mainly 
comes from the simplification of the spin vertex rather than 
the IPT approximation itself . We have checked that the Tn{U) 
curve becomes similar to the result of the QMC, if we include 
the (jL)n dependence in the spin vertex (neglecting off-diagonal 
elements of w„ and w^). However, in strong U region, the ac- 
tual Tn is also expected to be lower than that of DMFT with 
QMC, since DMFT completely neglects the spatial correla- 
tions. It is also curious to see how the nonlocal spin fluctu- 
ation affects the developed coherent quasiparticle state. This 
question will be answered in the next section. 

4. Influence of Spatial Correlation 

In this section we discuss the influence of the spatial cor- 
relation to the local self-energy based on eq. (2.8), in which 
local quantities are evaluated by DMFT (IPT). We have used 
64 X 64 discretized meshes in the 1st Brillouin zone and 2048 
Matsubara frequencies in the computation. 

Figure 8 shows the contour plot of the spectral intensity 
for U - % and T - 0.245 along the high-symmetry Unes of 
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Fig. 9. The density of states for C/ = 8. (a) temperature dependences, (b) 
comparison with the DMFT results. 



Fig. 10. The density of states for U = 14. (a) temperature dependences, (b) 
comparison with the DMFT results. 



the Brillouin zone. The overall structure of the spectral in- 
tensity is roughly the same as that obtained by eq. (2.4) with 
the DMFT self-energy (not shown). In contrast to the DMFT 
intensity, the pseudo-gap behaviors at w = appears in the 
coherent quasiparticle state due to the strong AF spin fluctu- 
ation (see, the X-M' line).^^^' Note that there exist the hot 
spots over the entire Fermi surface, since the magnetic Bril- 
louin zone coincides with the Fermi surface. The pseudo-gap 
phenomena due to the AF spin fluctuations are also discussed 
in the framework of the cluster DMFT.'^ 

To see how the pseudo gap develops with decrease of T, the 
total DOS, A(w) = -(1/A^) Zk ImG(k), for f/ = 8 is shown in 
Fig. 9. Figure 9(a) exhibits first the development of the co- 
herent quasiparticle peak at w = 0. As T approaches Tn =f 
0.241, A(w) yields the pseudo-gap behavior, and simultane- 
ously the upper and the lower Hubbard bands become more 
pronounced. The comparison with the DMFT DOS is shown 
in Fig. 9(b). For the temperature much above Tn, the inclu- 
sion of the spatial correlations affects little on A((l)), i.e., the 
spectral weight is slightly shifted to high-energy side. On the 
other hand, the quasiparticle peak is significantly suppressed 
near yielding the pseudo-gap behavior, and considerable 
amount of the weight in Hubbard bands is transferred to high- 
energy side. This spectral weight transfer brings the electronic 
state close to that in the AF ordered phase. "^^'^-^ 

The similar plot in the insulating case, U - 14 is shown in 
Fig. 10. In contrast to the results in the metallic case, the low- 
energy states near the insulating gap edge become more pro- 
nounced as T approaches Tn ==0.130. As similar to the metal- 



lic case, the weight of the Hubbard bands is shifted to high- 
energy side. However, its intensity is rather sharp as compared 
with that in the metallic case. In the present calculation we in- 
clude only the single exchange of the spin fluctuations. If we 
take account of multiple scattering from the spin fluctuations, 
the sharp incoherent peak would become much broader. Note 
also that the resultant spectral intensity near Tn has a close 
resemblance to that in the AF ordered phase. 

Finally, we discuss the T dependence of the specific heat, 
which is calculated by the numerical differentiation of the in- 
ternal energy. 



^=2j^£kG(k)e"''-°- +Ud. 



(4.19) 



Here the double occupancy, d is given by 

d^(lf +^J]Uk)G(k). (4.20) 

k 

The T dependence of the specific heat for f/ = 8 is shown 
in Fig. 11. The results by IPT are also shown for compar- 
ison. Owing to the short-range AF order, the specific heat 
is strongly enhanced toward Tn. As shown in the inset, the 
double occupancy rapidly decreases toward the AF instabil- 
ity^^ in contrast to the increasing behavior in IPT. It should be 
noted that the strong enhancement in the specific heat near Tn 
mainly comes from the reduction of d rather than that of the 
kinetic energy. In other words, the AF correlation enhances 
the tendency of the Mott localization. The QMC simulation 
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Fig. 1 1 . The T dependences of the specific heat for t/ = 8. The inset shows 
the T dependence of the double occupancy. The results of DMFT is shown 
for comparison. 



using the 6 X 6 square lattice shows the T dependence of the 
specific heat similar to that of the DMPT."*"* This implies the 
importance of the long-range critical fluctuations. 

5. Summary and Discussion 

We have proposed the formulation to take into account the 
correction of the spatial correlations to the local self-energy 
obtained by DMFT. It is demonstrated that the one-particle 
spectral intensity of the 2D Hubbard model at half-filling 
shows the pseudo-gap behavior in the central coherent quasi- 
particle peak due to the critical AF fluctuation. It is important 
to note that the pseudo gap driven by the critical AF spin fluc- 
tuation is always accompanied by the enhancement ofl/TiT 
in the NMR relaxation rate. Thus, the spin-gap behavior-'"* in 
l/TiT of the hole-doped high-Je cuprates in the underdoped 
region requires alternative explanations such as the eff'ect of 
critical superconducting fluctuations.'*^'"*^ The specific heat is 
considerably enhanced by the AF short-range order, which is 
accompanied by the reduction of the double occupancy, indi- 
cating the enhanced tendency of the Mott localization. 

Here we should note a deficiency of the approximate lo- 
cal dynamical susceptibiUty based on IPT. The spectral in- 
tensity of the approximate dynamical spin susceptibility be- 
comes slightly negative in high-energy region. This is because 
the use of the energy-independent charge vertex in eq. (3.17). 
Note that the charge vertex should become the bare U in 
the high-energy hmit. This deficiency can be removed by 
the use of a charge vertex with proper energy dependence, 
which is evaluated by QMC for instance. The QMC simula- 
tion with the Hirsch-Fye algorithm is of course more accurate 
than the present method, however, the computational cost is 
much more expensive, especially in lower temperatures and 
the strong U region. It should be emphasized that the defi- 
ciency of the present method is not so serious in the improved 
self-energy, since we consider only the nonlocal correction of 
the susceptibihties, h-Xaiq) = Xaiq) - Xai^m), where Axaiq) 
vanishes in high-energy region as it should. In other words, 
the essential local dynamics is sufficiently taken into account 
in the DMFT self-energy. 



Throughout this work, we have updated the self-energy 
only once. Therefore, the spin fluctuation in our self-energy 
contains no feedback from the formation of the pseudo gap. If 
such a feedback is taken into account, the critical AF fluctua- 
tion as weU as the pseudo-gap behavior would sUghtly be sup- 
pressed. In order to consider the feedback effect, one possible 
way is to solve the Dyson's equation in a self-consistent man- 
ner, keeping the local vertices and the DMFT self-energy un- 
changed. This kind of self-consistent treatment includes mul- 
tiple scattering from spin fluctuations, which would make the 
incoherent part of the single-particle spectrum much broader. 

A further improvement of the present formulation should 
consider the local moment sum-rule (the fluctuation- 
dissipation theorem) for the spin susceptibihty,*" which is 
slightly violated in the present formulation. This sum rule be- 
comes significantly important to discuss the anomalous tem- 
perature dependences in the vicinity of the quantum critical 
point."*^^^ In view of the success of the self-consistent renor- 
malization theory,"*^- an introduction of a constant shift in 
the static component of the spin irreducible vertex may work 
out to satisfy the local-moment sum-rule. 

On the physics of the 2D Hubbard model in the context 
of the high-Tc cuprates, it is important to discuss the effect 
of an asymmetry of the tight-binding dispersion with longer- 
range hoppings, and of electron/hole doping.^* In region of 
finite dopings, it is also curious to see the superconducting 
instability based on the formula given in Appendix B. They 
are left for the future work. 
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Appendix A: Iterated Perturbation Theory 

In this Appendix, we summarize the formulation of the it- 
erated perturbation theory (IPT) as a DMFT solver. Motivated 
by the success of the second-order perturbation theory in the 
symmetric impurity Anderson model,^' the local self-energy 
in IPT is written in the interpolation form,^^'^'^^* 

A2:(2)(a;„) 



l-B5:P)(a>„)' 
where the second-order self-energy is expressed as, 



(A-1) 



(A-2) 



wifli A'o(em) = -T Zco„ @i(^n)@ipJn + u)- &ipJn) IS file cavity 
Green's function or the "bare" Green's function in the ficti- 
tious impurity Anderson model, which contains aU of the dy- 
namical information of the other sites of the lattice. 

The coefficients, A and B, are determined so as to satisfy 
the high-energy and the atomic Umits.^^^^ Their expressions 
are given by 



B = A 



n(2 - n) 
noil-no)' 
jlf-fi + U{1 - ri) 



(A-3a) 



(A-3b) 



f/2no(2 - no) ' 
where the electron density n = 2TY,o> G^(w„)e'""°* and 
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riQ = 2T 2(j„ ^(w„)e"^"*'+ . Here the fictitious chemical poten- 
tial of the mapped Anderson model, jlf, is introduced, which 
is determined by the condition in the low-energy limit.^^ 
Namely the Luttinger sum-rule is enforced to determine jlf 
at r = 0.^'''^^ For finite temperatures, a natural extension of 
the Luttinger sum-rule is the following identity,'' 



dia>n 



= 0, 



(A4) 



which is evaluated via the analytic continuation with the Pade 
approximant. 

Using the local self-energy, (A-1), we obtain the local 
Green's function as 

G\oJn) = ^.Ycik) ^ r /"ff, , . (A-5) 



where po(e) is the bare DOS of the conduction electron. The 
cavity Green's function is calculated by 



Then, the updated Q{iOn) is used to calculate new self-energy 
until a convergence within the given accuracy is reached. Note 
that in the symmetric case,^^ we have simply ftf - jl = and 
„ = „o = 1 (A = 1 and B = 0), and tHoJn) = 2:P)(w„). 

There is one technical comment in order to obtain an in- 
sulating solution. The insulating solution has the form in the 
limit of |w„| ^ 0,22-23 



(<^n) ~ 1- 2inc» 

ico„ 

b ItOn 



b 



(A-7a) 



(A-7b) 



(A-7c) 



The last expression with the leading term i(jL>nlb is obtained by 
the exact cancellation of the singular terms, ajiojn in fy{<jJn) 
and G~^{tL)„). In practice, however, numerical errors often 
leave the singular term a/Zw,, in ^(w,,), which masks the lead- 
ing term iojn/b. This results in the convergence to the metallic 
solution instead of the insulating one. To avoid miscancella- 
tion, eq. (A-5) is modified to 



z„ 



I' 



1+ de 



(A-8) 



withz„ = ia)„+jl-'l}^ia)„). Then, the factor l/z„ is responsible 
for the exact cancellation with f^{(o„). 

Appendix B : Determination of Superconducting Transi- 
tion Temperature 

Within the approximation proposed in this paper, the super- 
conducting transition temperature, Tc, can also be determined. 
In this appendix, we describe how to determine Tc for the sin- 
glet and the triplet channels. 

In general, Tc is determined by the divergence of the 
static and homogeneous component of the renormalized ver- 
tex function, r*(fc, fc'; 0)."* T* obeys the BS equation in the 
particle-particle (PP) channel, 

r;(A:,r;0) = rf(A:,A:';0) 







,-k* 


k'- 




P-, 






k. 




>-k 


k' 



k"+k'-k k"'-Hk'-k 

'u" 



-k' 



local 
vertex 



Fig. B-l. The Feynman diagram of the particle-particle irreducible vertex 
expressed by the spin and the charge fluctuations. 



-z 

k" 



(A:, k"; 0)\Gik")f-r^ik",k'; 0), 



(Bl) 



where F^^ (77 = sin, tri) is the irreducible vertex in the singlet 
and the triplet PP channels. In the vicinity of the transition, 
we can neglect the first term to yield 

= -AY,T^''ik,k'-MGik')f-Laik'), (B-2) 

k' 

where the anomalous self-energy is introduced as £a(A:) = 
(A-6) ~ ^n^^' Then, Tc is determined such that the 



maximum eigenvalue A of the above eigenvalue equation 
becomes unity. In practice, the alternative expression of an 
eigenvalue problem for the anomalous propagator is useful, 

,p,(k) = -A\G(kf Yj rf^^' ^)<Pn(k'\ (B-3) 

k' 

since (f>ri(k) decays faster than 1 /\co„[^ in the high-energy limit. 

To solve the eigenvalue equation, it is necessary to give 
the expUcit form of the irreducible vertices. In general, they 
are given by the functional differentiation of the Luttinger- 
Ward functional with respect to the anomalous Green's func- 
tion,^' ^' 2^ and they can be expressed in terms of the irreducible 
vertices in the particle-hole channel. When the superconduct- 
ing transition is triggered by the spin and/or the charge fluc- 
tuations, the dominant contribution to the irreducible vertices 
is given by the exchange of the single collective modes (so- 
called the Maki- Thompson process).^' ^'^^ Neglecting con- 
tributions other than the Maki-Thompson process, the irre- 
ducible vertices can be expressed within the present approxi- 
mation as 



-I 



-3r^^(co'„-co„)Axsp(k'-k) +FLja>„,w;;0), (B4a) 

J even 

'yii,(co'„-aJn)AXcb(k'-k) 



r^^(k,k';0) = - 



odd 



(B4b) 



where the subscript even (odd) indicates the even (odd) part 
with respect to ^ — > -k. Here, F^j^^ is associated with the local 
anomalous self-energy, Ej;'(fa;„). Its explicit form is not nec- 
essary in considering the anisotropic pairings, since F'^ does 
not affect the transition temperature for anisotropic pairings. 
The Feynman diagram of the PP irreducible vertex is shown 
in Fig. B l. 
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